Physiological Genomics

Maternal hyperglycemia can affect offspring development and metabolism, but how its effects on skeletal muscle change during postnatal life remains unclear. We investigated skeletal muscle phenotypes in male offspring of streptozotocin-induced hyperglycemic mice at 4 and 8 weeks of age, representing early post-weaning and later postnatal stages, respectively. Maternal hyperglycemia resulted in lower offspring body and skeletal muscle weights at both time points. Myofiber cross-sectional area was reduced without a reduction in total myofiber number, indicating impaired myofiber growth rather than myofiber loss. At 4 weeks, skeletal muscle showed type I/IIA hybrid myofibers and increased expression of fetal, myogenic, and proliferation-associated genes. This immature profile was not maintained at 8 weeks; instead, the transcriptional profile shifted toward reduced mitochondrial and oxidative metabolism. DNA methylation changes at 4 weeks showed limited correspondence with transcriptional alterations. Taken together, maternal hyperglycemia impaired postnatal myofiber growth at both time points, with distinct molecular profiles at 4 and 8 weeks.

ABSTRACT Objective: Skeletal muscle growth rate differs substantially among diverse pig breeds. In the present study, RNA-seq was performed to identify differentially expressed genes (DEGs) across pig breeds with distinct growth rates. CCND2 showed obvious differential expression among pig breeds. Although CCND2 is essential for cell cycle control, its biological function in skeletal muscle differentiation has rarely been reported. This work aims to identify DEGs in pig breeds with different growth rates and explore CCND2 role in porcine skeletal muscle development. Methods: RNA-seq was performed to identify differentially expressed genes (DEGs) across pig breeds with distinct growth rates. Weighted gene co-expression network analysis (WGCNA) was used to identify gene modules potentially associated with skeletal muscle growth and development. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis was performed to validate the role of CCND2 in skeletal muscle differentiation. Results We identified 5,416, 4,846, and 2,777 DEGs in the Duroc pig (DU) vs. Jiangquan black pig (JQ), DU vs. Laiwu pigs (LW), and JQ vs. LW comparison groups, respectively. WGCNA identified two gene modules (lightgreen and darkmagenta) potentially associated with skeletal muscle growth and development. CCND2 exhibited differential expression in comparisons between DU vs JQ and DU vs LW, with higher expression observed in Duroc pigs. qRT-PCR analysis revealed that overexpression of CCND2 at day 6 (D6) of cell differentiation significantly increased the expression levels of CCND2 gene and markers of myogenic differentiation (p < 0.05). Conversely, knockdown of CCND2 significantly decreased the expression levels of CCND2 gene and myogenic differentiation markers (p < 0.05). Conclusions This study identified breed-specific DEGs and key gene modules in porcine skeletal muscle and verified that CCND2 promotes myogenic differentiation. These findings provide a genetic theoretical basis for further investigations into skeletal muscle growth and development. Keywords pig, skeletal muscle, transcriptome, WGCNA, CCND2

Introduction: Adipose tissue fibrosis has been associated with metabolic dysfunction-associated steatotic liver disease (MASLD) and suboptimal weight loss after bariatric surgery. The objective was to characterize changes in adipose tissue fibrosis and extracellular matrix (ECM) components in adipose tissues from two abdominal fat compartments, subcutaneous and omental, in response to bariatric surgery; and determine the effect of weight loss on the association between adipose tissue fibrosis and MASLD feature changes with weight loss. Methods: Liver and adipose tissue biopsy samples were obtained before and after weight loss from patients who underwent two-stage bariatric surgery. Adipose tissue pericellular fibrosis was evaluated histologically with picrosirius red staining for collagen quantification and normalized to adipocyte count per microscope field. Transcript levels of ECM remodeling and inflammatory genes were measured by RNA sequencing in subcutaneous adipose tissue pre- and post-surgery. Results: At baseline, adipose tissue pericellular fibrosis did not associate with MASLD features, metabolic dysfunction-associated steatohepatitis or hepatic fibrosis in patients with severe obesity. Changes in liver histology after surgery were unrelated to changes in adipose tissue fibrosis. Transcripts for most genes involved in ECM remodeling and inflammation were downregulated in subcutaneous fat one year after bariatric surgery. Omental, but not subcutaneous, adipose tissue fibrosis decreased after surgery when adjusted for cell number (mean {+/-} SD: 0.41 {+/-} 0.12 vs 0.32 {+/-} 0.08 {micro}m2/cell, p<0.005). Higher subcutaneous adipose tissue fibrosis correlated with reduced weight loss (estimate {+/-} SE: -2.04 {+/-} 0.99, p<0.05) and BMI loss (-0.73 {+/-} 0.35, p<0.05), while increased omental adipose tissue fibrosis after weight loss was associated with worsening hepatic fibrosis (0.31 {+/-} 0.14, p<0.05). Conclusion: In individuals with severe obesity undergoing bariatric surgery, adipose tissue fibrosis is poorly correlated with MASLD features. Notably, pericellular fibrosis in subcutaneous abdominal, but not omental, adipose tissue was associated with a reduced weight-loss response to surgery.

Sepsis mimics resemble sepsis clinically but may require different treatments. We evaluated whether inflammatory protein profiles could distinguish these conditions in a retrospective case-control study using 4 prospective Japanese cohorts. The study included 34 adults with sepsis and 46 with sepsis mimics. We measured 92 proteins using Olink, assessed between-group differences with false discovery rate (FDR) adjustment, and evaluated 3 prediction models using repeated nested cross-validation. No protein met the FDR threshold; the lowest FDR-adjusted P value was 0.345. The primary elastic-net model had an area under the receiver operating characteristic curve of 0.442 (95% confidence interval, 0.315-0.574), based on cross-validation predictions averaged for each participant. Random forest and gradient boosting also showed no clear evidence of discrimination above chance. Confidence intervals included chance performance in both the all-serum and quality-control-Pass sensitivity analyses. In this cohort, a single measurement of an inflammation-focused protein panel showed limited discrimination between sepsis and sepsis mimics. Integration with clinical information and longitudinal measurements warrants further investigation.

Postmenopausal women are vulnerable to metabolic dysfunction. Compelling evidence suggests this is because they lack the protective effects of estrogens. We have shown that circulating estrogens regulate daily eating and activity rhythms in female mice and protect them from obesity and metabolic dysfunction. However, few studies have investigated whether postmenopausal women have disrupted eating and sleep-activity rhythms that could contribute to metabolic dysfunction. The purpose of this cross-sectional study was to investigate the relationship between sleep and eating rhythms and metabolic risk in postmenopausal women with overweight or obesity. Seven days of sleep and food intake behavior data were collected from postmenopausal women with overweight or obesity. Sleep timing was assessed using wrist-worn actigraphy and sleep logs. First and last calorie intake times were collected daily using text messaging. Visceral adipose tissue (DXA), insulin resistance (Matsuda index), blood pressure, and triglycerides were analyzed as markers of metabolic risk. In 122 postmenopausal women (mean {+/-} SD: 58.1 {+/-} 3.6 yr, 33.2 {+/-} 5.2 kg/m2 BMI), we found that later sleep onset was associated with greater visceral adipose tissue (p=0.0007), but not insulin resistance, blood pressure, or triglycerides. Next, we used a random forest model to quantify the relative contributions of sleep and meal timing variables to predict visceral adipose tissue. The top predictors of visceral adipose tissue were time of sleep onset, sleep duration, and time of last caloric intake. These data suggest that interventions that advance sleep onset may reduce metabolic risk in postmenopausal women.

BACKGROUND Twin pregnancies have a two- to three-fold higher risk of preeclampsia (PE), but early prediction remains poor because current biomarkers and prediction models are largely derived from singleton pregnancies. Defining the biological changes that precede PE in twins may facilitate the development of more effective risk assessment approaches. OBJECTIVE We aimed to characterize maternal plasma protein changes preceding PE in twin pregnancies and to explore whether early circulating proteins could serve as improved risk stratification for twin PE. STUDY DESIGN We performed a nested case-control study within a prospective twin pregnancy cohort at Peking University Third Hospital. A total of 84 women with twin pregnancies were included, comprising 40 pregnancies complicated by PE and 44 normotensive controls. Longitudinal maternal plasma samples were collected in the first, second, and third trimesters. In the discovery cohort, maternal plasma proteins were profiled by mass spectrometry to identify trimester-specific protein changes associated with PE. Candidate proteins were then validated in an independent validation set. We further assessed the functional relevance of selected proteins in trophoblast and endothelial cell models and constructed prediction models integrating biomarkers with maternal clinical characteristics. RESULTS Twin pregnancies complicated by PE showed distinct trimester-dependent protein alterations, with first- and second-trimester changes mainly involving immune regulation and glycolipid metabolism, and later changes enriched in cardiovascular-related pathways. Several proteins demonstrated early discriminatory patterns, including increased annexin A3, chitotriosidase-1, and ferritin light chain (FTL), together with decreased annexin A2 in the first trimester, whereas the conventional singleton markers (sFlt-1 and PAPP-A) did not differ. ELISA validation of these candidate biomarkers yielded an AUC of 0.93 for twin preeclampsia prediction; integration with standard maternal clinical factors further improved discriminatory performance to an AUC of 0.97. Accordingly, a multivariable model combining annexin A2, annexin A3, chitotriosidase-1, FTL, and maternal clinical factors achieved strong internal performance for early discrimination of PE in twin pregnancies (AUC, 0.97; 95% CI, 0.93-1.00). Moreover, annexin A2 knockdown impaired trophoblast migration and invasion and reduced inflammatory cytokine secretion in endothelial cells, supporting a role in placental development and maternal vascular adaptation. CONCLUSION Twin pregnancies that later develop PE exhibit early and progressive maternal plasma protein alterations that can be detected before clinical disease onset. These findings identify candidate biomarkers and biologically relevant pathways that may help explain why PE in twins differs from singleton disease and may support future twin-specific risk stratification strategies. These results require external validation before clinical application.

Given unprecedented environmental change, Pacific salmonids are becoming increasingly threatened; however, genomic technologies provide an opportunity to assess Pacific salmonid resilience to on-going environmental change. Also, given increasing conservation and animal welfare concerns, non-lethal tissue sampling is emerging for assessing the physiological status of wild species. Relatively few studies, however, have compared the transcriptomic response of fish sampled non-lethally versus lethally for fish health and stress quantification. Here, we conducted a comparative gene transcriptional profiling study using 28 qPCR assays on the OpenArray platform as an osmoregulatory stress transcriptional profiling chip. We sampled both non-lethal (gill biopsies, skin swabs, and water [eRNA]) and lethal (whole gill and liver) tissue samples from Chinook salmon (Oncorhynchus tshawytscha) smolts exposed for 0 h, 24 h, and 72 h to freshwater (0 PSU), brackish water (15 PSU) or seawater (30 PSU). We observed that most assays showed higher qPCR amplification success in both non-lethal gill biopsies and lethal gill tissues, compared to other tissues. We found more significant differences in transcript abundance in traditional lethal tissues compared to non-lethal tissues. With increasing duration of salinity exposure, gene transcript levels were typically upregulated in seawater samples compared to freshwater samples for both non-lethal and lethal tissues. The targeted genes for both non-lethal and lethal tissues were significantly and positively correlated across time points. Our findings support the use of non-lethal tissue sampling to improve Chinook salmon hatchery rearing, inform hatchery release timing, and facilitate future field application for the conservation of wild stocks.

Obesity and metabolic diseases such as diabetes reflect long-term interactions between diet, lifestyle, and genetics, variables that are hard to control in humans but can be systematically varied in mice, though behavioral variables like long-term exercise remain harder to model. Having previously phenotyped ~50 BXD strains, we identified BXD40, BXD43, and BXD100 as predisposed to voluntary wheel-running (~6 km/day at 24 weeks) with variable susceptibility to diet-induced obesity and diabetes. We followed these strains from 8-32 weeks on a 60% kcal/fat high-fat (HF) diet, with or without access to running wheels. Wheel usage declined from ~9 km/day in the first two weeks to a steady ~5 km/day thereafter yet still improved oral glucose tolerance by 35% over sedentary controls by 24 weeks. Despite equivalent body weight, runners ate 20% more relative to body weight, indicating the glucose benefit was independent of caloric intake. We next asked whether the gut microbiome contributes to this benefit. Long-term exercise produced only modest, cohort-dependent shifts in cecal microbial composition, though functional profiling revealed exercise-associated differences in carbohydrate- and glucose-metabolism pathways. Reciprocal bedding exchange between sedentary and running BXD43 cages reduced microbial differentiation in a taxon-specific manner and coincided with marked attenuation of the glucose phenotype, despite unchanged running behavior. Together, these findings validate BXD40, BXD43, and BXD100, alongside the canonical C57BL/6, as a model for studying diet, exercise, and genetics in metabolic disease, and implicate gut microbial taxa in shaping exercise's metabolic effects during HF feeding.

Background: Dietary sodium intake is a modifiable determinant of blood pressure, and impaired renal sodium handling is implicated in the pathogenesis of hypertension. Whether dietary sodium intake modulates renal tubular sodium transporter abundance among normotensive adults at risk for hypertension is unknown. Methods: Forty-six normotensive adults with normal kidney function at risk of developing hypertension were prospectively enrolled to complete a very low sodium (LS) diet followed by a high sodium (HS) diet. Urinary extracellular vesicle (uEV) proteins were quantified from a 24-hour urine collection at the end of each diet. We compared normalized uEV protein concentrations between diets and tested whether diet modified uEV protein associations with serum aldosterone, urinary aldosterone, and serum cortisol-cortisone ratio under either diet using Wilcoxon rank sums and linear mixed effects models, respectively. Results: Participants achieved median [IQR] urine sodium 13.2 [7.7-19.9] mEq/24h on the LS diet and 244.0 [197.1-297.0] mEq/24h on the HS diet. LS intake induced greater renin activity, greater plasma and urinary aldosterone, and lower urinary cortisol, when compared to HS intake. In parallel, NCC, ENaC, SGLT2, OXSR1, and Pendrin uEV concentrations were significantly higher during LS intake than HS intake (p<0.05 for all), whereas TRPV5 and AQP2 uEV concentrations were lower during LS intake than HS intake (p<0.001 for both). Conclusions: Controlled dietary sodium modulation induced nephron-segment specific changes in uEV protein abundance. These findings demonstrate the ability of uEVs to serve as a non-invasive tool to probe in vivo renal sodium, water, and electrolyte handling.

Synchrony between time-of-day and food intake is critical for maintaining homeostasis while desynchrony between these cues caused by repeated bouts of mistimed eating (ME) is profoundly disruptive to organismal health. However, it is unknown whether the same is true for a single bout of ME. While it is generally accepted that mammals recover quickly from acute alterations in feeding schedule and that chrono-disruption requires chronic, repeated changes, there have been no dedicated investigations of acute ME to suggest otherwise. Therefore, we exposed female and male C57Bl6/J mice to a single bout of ME or timed eating (TE) and continuously measured metabolic and behavioral responses for five days. A single bout of ME was sufficient to disrupt the daily rhythms of metabolic and behavioral outcomes for several days. Energy expenditure and wheel running activity were disrupted and light:dark variation in food intake and macronutrient utilization were altered for several days. Finally, ME changed the amount of time spent resting in for several days post-ME. The TE paradigm disrupted macronutrient utilization for two days but otherwise had no effect on metabolism or behavior. Our results demonstrate that one instance of ME alone has lasting impacts on daily rhythmicity and is an understudied source of circadian disruption.

Objective: To test the hypothesis that selected SIRS and SOFA criteria measured at ICU admission have divergent associations with (1) the development of new culture-confirmed infections and (2) early empiric broad-spectrum antibiotic receipt in critically ill patients with cirrhosis. Design: Retrospective, observational study using electronic health record (EHR) data. Setting: Both community and academic ICUs within a five-hospital health system between July 2017 and December 2025. Subjects: Adult patients >18 years of age diagnosed with cirrhosis admitted to the ICU within 96 hours of hospital admission. Patients transferred from outside facilities and those with positive cultures collected >24 hours before ICU admission were excluded. Interventions: None. Measurements: Multivariable logistic regression tested the association of the most aberrant values in nine SIRS/SOFA vital sign and laboratory parameters measured during the 24 hours prior to ICU admission with (1) the development a newly positive culture collected during the 24 hours before through the 96 hours after ICU admission and (2) receipt of broad-spectrum antibiotics during the 24 hours before through the 24 hours after ICU admission. Main Results: The presence of hyperthermia, hypothermia, hypoxia, or hypotension at ICU admission was associated with both culture positivity and antibiotic receipt. The presence of leukocytosis was associated with antibiotic receipt but not culture positivity. Among patients in the 21.6% of encounters with newly positive cultures, normothermia and normal WBC were common (respectively present in 97.2% and 77.5% of encounters) and associated with antibiotic omission at ICU admission despite the subsequent development culture positivity. Conclusions: Traditional physiologic markers of infection are often absent in critically ill patients with cirrhosis who develop new culture-confirmed infections at ICU admission, possibly contributing to delayed antibiotics. There may be opportunities to improve antibiotic delivery and stewardship at ICU admission in this high-risk population.

Type 2 diabetes mellitus (T2DM) is a chronic disorder associated with obesity and hyperglycemia. T2DM is linked to significant clinical complications such as delayed wound healing and diabetic foot ulcers. Tumor necrosis factor-stimulated gene 6 (TSG-6) is a protein with multiple functions including regulation of cytokine binding to glycosaminoglycans in extracellular matrices and on vascular endothelial walls. Loss of TSG-6 in knockout mice increases skin inflammation and delays wound healing. Here, in both a genetic model (db/db mice) and a dietary model that simulates human T2DM (mice treated with a high-fat diet [HFD] and streptozotocin [STZ]), we observed sustained hyperglycemia and a reproducible delay in wound closure. TSG-6 expression was markedly reduced in diabetic skin in these models, both before injury and post wounding. HFD+STZ diabetic mouse wounds showed elevated levels of pro-inflammatory cytokines (IL-6, TNF-, and MIP-1) and lower levels of the anti-inflammatory cytokine IL-10, relative to normoglycemic mice. Re-introduction of TSG-6 protein into the wounds normalized healing rates, corrected abnormal cytokine expression levels, and increased macrophage accumulation in dermal and adipose compartments. Thus, T2DM is associated with reduced levels of TSG-6 in diabetic skin, which contributes to the pro-inflammatory milieu and the delay in wound closure.

Environmental exposures can induce epigenetic modifications that persist across generations, potentially contributing to the transmission of environmentally-induced phenotypes. However, the extent to which such persistence of molecular changes are maintained independently of genetic mechanisms remains unclear. In this study, we investigated the evolution and genetic determinism of DNA methylation (DNAm) variation across three successive generations of Japanese quails (G0 to G2) belonging to two epilines, defined according to whether their female ancestor had received genistein supplementation (epi+) or not (epi-). Using reduced representation bisulfite sequencing (RRBS), we characterised CpG methylation patterns and identified differentially methylated cytosines (DMCs) and regions (DMRs) between control and genistein-supplemented epilines. Out of 112,745 CpG sites analysed, no DMCs or DMRs were detected between epilines in the first generation following exposure (G0), whereas 621 DMCs (36 DMRs) were identified in G1, and 1,381 DMCs (101 DMRs) in G2. Despite this progressive increase of differential methylation sites, mean differences in methylation rate between epilines were globally small. Similarly, genome-wide genetic differentiation between epilines was very limited. Overall, SNP-based heritability of DNA methylation was 0.19 indicating that genetic variation accounted for only a small part of the global variation of DNAm levels. By contrast, DNAm at differentially methylated CpG sites was extensively under genetic control, with a mean SNP-based heritability > 0.5, and meQTL analysis further identified significant associations between SNPs and methylation levels. Together, these results indicate that local DNAm differentiation in quail is under substantial genetic regulation, while the limited genetic differentiation between groups suggests that genetic variation alone can not fully explain the progressive accumulation of methylation differences following an ancestral exposure to genistein. These results support DNAm as a plausible molecular candidate mechanism in the multigenerational response to environmental exposure and highlight the complex interplay between genetic and epigenetic regulations in the inheritance of phenotypes.

Enhanced excitability of output neurons in the olfactory bulb (OB) has been linked to regulation of energy homeostasis. The mechanism by which the olfactory system can drive metabolic state while processing odor information remains elusive yet must depend upon integration and relay to the rest of the brain. Herein, we screened brain-wide neuronal ensembles that preferentially engaged during odor stimulation (referred to as ODORTRAP) utilizing FOSTRAP activity mapping. Mice were exposed to the odor isopropyl tiglate (IPT), a preferred ligand for the olfactory receptor 160, previously shown to be vulnerable to diet-induced obesity. Odor activation resulted in elevated labelling in the olfactory bulb (ODORTRAP) (p < 0.001), as well as in higher olfactory cortices including anterior olfactory nucleus, piriform cortex, and cortical amygdala (p < 0.05). Interestingly, IPT differentially recruited hypothalamic feeding areas, namely the arcuate (ARC), dorsomedial (DMH) and lateral (LH) hypothalamus (p < 0.001), while the ventromedial hypothalamic nucleus and the suprachiasmatic nucleus exhibited no change in ODORTRAP activity (p > 0.05). In mice that were maintained on moderately high-fat (MHF) diets, we discovered that obesity re-models the activation pattern across the brain. Obese mice exhibited loss of ODORTRAP activity in the ARC (p < 0.0001), whereas in the DMH and LH increased baseline activity that obscured any odor-linked neural activity. While IPT activated the nucleus accumbens and medial septum in control mice (p < 0.01), obese mice demonstrated no changes in ODORTRAP activity (p > 0.05). Our study identifies key brain regions linking chemosensation and metabolic state, specifically demonstrating alteration in activity in hypothalamic domains upon obesogenic diet.

Objectives: The bile acid nuclear receptor Farnesoid X Receptor (FXR, NR1H4) is a major regulator of metabolism and energy homeostasis in peripheral organs. It modulates bile acid, glucose, and lipid metabolism, as well as fat mass and body weight. However, FXR is also expressed in the brain, particularly in the hypothalamus, a key center for the regulation of energy homeostasis. Although one study has demonstrated a role for brain FXR activation in energy balance, its specific hypothalamic role is still unknown. Here, we examined the role of FXR in the mediobasal hypothalamus in the regulation of energy balance. Methods: We used a genetic approach combined with metabolic phenotyping to determine the effect of FXR invalidation in the mediobasal hypothalamus on metabolic parameters involved in the central regulation of energy homeostasis. Results: Our results demonstrate that hypothalamic FXR deficiency induces a positive energy balance, resulting in a reduction in energy expenditure due to alterations in glucose metabolism accompanied by structural changes in white adipose tissues. Conclusion: This study uncovers a previously unrecognized role for hypothalamic FXR in the central homeostatic control of energy balance, providing new insights into its contribution to peripheral glucose metabolism and adipose tissue structural remodeling.

Genetic loci influencing multiple complex traits often remain difficult to interpret because the underlying functional variants, target genes and tissue contexts are unresolved. Here, we investigated the regulatory basis of a shared signal for average daily gain (ADG) and backfat thickness (BF) in Yorkshire pigs. By integrating GWAS in 7,972 pigs, targeted capture sequencing, LD-based fine-mapping, porcine epigenomic annotations, three-dimensional chromatin interaction maps and functional assays, we narrowed the shared SSC1 signal to a 1.59-Mb interval and functionally prioritized rs334769214 (T>G) as the leading candidate variant. The T>G substitution produced marked allele-specific differences in enhancer activity and nuclear-protein binding, while transcription-factor co-transfection assays identified FOXL1-responsive, allele-dependent regulatory activity. Chromatin interaction and CRISPRa analyses supported a regulatory connection between the rs334769214-containing enhancer and PMAIP1, and PMAIP1 overexpression promoted apoptosis in porcine intestinal epithelial cells. In the capture-sequenced subset, the G allele was associated with higher ADG and BF, consistent with a relative loss-of-activation model. Together, these findings delineate a candidate FOXL1-enhancer-PMAIP1 regulatory pathway and demonstrate how multi-omics fine-mapping can advance a shared quantitative-trait locus toward a testable molecular mechanism underlying variation in growth and fatness.

Stem cell-derived (SC) islets hold promise as a transformative therapy for type 1 diabetes (T1D). To enhance graft uniformity, SC-islet organoids are often dissociated and reaggregated into standardized sizes. However, how endocrine cell composition varies with islet size in native human islets remains largely unknown. Defining this relationship is essential for establishing physiological benchmarks for islet engineering for transplantation. In this study, we examined the relationships among islet size, endocrine cell composition, and spatial distribution throughout whole pancreata from 10 age-matched, non-diabetic male subjects. Insulin-, glucagon-, and somatostatin-positive endocrine cells were identified as single-, double-, and triple-hormone populations ranging from non-islet single endocrine cells and small clusters (niSSCs) to large islets. Approximately 80% of the total endocrine area was contained in islets positive for all three hormones. Islet size strongly correlated with endocrine composition. Specifically, larger islets had increased alpha-cell fractions and reduced beta-cell fractions, while delta-cell fractions were consistent across sizes. Spatial analysis revealed that niSSCs, mono-, and tri-hormone islets frequently exhibit a clustered distribution. Additionally, tri-hormone islets were positioned near the core of the pancreas, whereas mono-somatostatin islets were often localized in the periphery. These findings were consistent across all three pancreatic regions (head, body and tail).

The transcriptomic interaction between altered gravity and innate immunity remains largely unexplored. In this study, we established a novel dual-stress model combining 3g hypergravity with lipopolysaccharide (LPS) immune challenge in 5 days post-fertilization (dpf) zebrafish (Danio rerio) larvae. RNA sequencing revealed that exposure to hypergravity alone induced no differentially expressed genes (DEGs), indicating a remarkable transcriptional resilience. In contrast, the addition of an immune challenge under hypergravity conditions triggered a substantial transcriptomic response involving 2,706 DEGs. Functional enrichment analyses identified significant alterations in pathways related to immune function, development, and circadian regulation. Targeted qPCR analyses further validated the expression patterns of selected genes. These findings suggest that while zebrafish larvae can effectively tolerate hypergravity alone, the combination of gravitational and immune stressors elicits profound molecular responses. This study provides new insights into the interplay between altered gravity and innate immunity using an animal model, contributing to future strategies aimed at safeguarding astronaut health and supporting space-based aquaculture systems.

Genetic variation impacts biological response to chemical and environmental exposures, contributing to differences in resistance or susceptibility. Forward genetic screens of genetically diverse cell populations can be used to identify the precise genetic variants that drive response variation. Cell lines from laboratory mouse genetic reference populations like the Diversity Outbred (DO) population are well powered for genetic mapping. Diverse stem cell panels are especially beneficial for studying early developmental effects of chemical exposures like triphenyl phosphate (TPHP), an organophosphate flame retardant linked to adverse developmental effects, including altered cell cycle in stem cells. Self-renewal and pluripotency are defining properties of stem cells during development; therefore, it is important to understand how chemicals like TPHP may alter these critical characteristics. To better understand the influence of genetic variation on TPHP exposure response, we utilized a genetically diverse panel of DO induced pluripotent stem cells (DO iPSCs). Our analyses characterize interline variation in gene expression through differential gene expression analysis and expression quantitative trait locus (eQTL) mapping. We identified genomic loci that contribute to variation in gene expression following TPHP exposure, including a regulatory hotspot on chromosome 15. Gene set enrichment highlighted several pathways affected by TPHP including lipid metabolism, steroid hormone signaling, and cell cycle. Our study shows that genetic variation modulates the effects of TPHP on gene expression and cell cycle in stem cells. Additionally, our work demonstrates that genetically diverse cell panels like the DO iPSC resource offer a tractable and scalable approach for identifying the genetic determinants of chemical toxicity.

The turquoise killifish (Nothobranchius furzeri) is an emerging vertebrate model for development and ageing due to its exceptionally rapid maturation and short lifespan. However, how microRNA (miRNA) regulation coordinates gene expression across its compressed developmental trajectory remains poorly understood. Here, we integrated miRNA and mRNA expression from the same unexposed animals across three developmental stages: newly hatched fry (4 weeks total age), juveniles (~8 weeks total age), and aged adults (~8 months of age). Developmental stage was the dominant source of variation in both datasets, with miRNA and mRNA profiles independently resolving a consistent progression from newly hatched fry to juveniles and adults. Transcriptomic changes reflected ongoing tissue maturation, including increased expression of crystallin genes associated with lens development, while more than 2,000 genes were differentially expressed during the juvenile transition. Integration of 398 developmentally dynamic miRNAs with predicted mRNA targets identified specific miRNAs whose developmental changes were significantly enriched for opposing changes in their target genes. These included conserved developmental regulators from the miR-200/141, miR-29, let-7, and miR-34 families, which increased toward adulthood as their predicted target cohorts decreased. Stage-controlled analyses further identified miR-34a, miR-29a, miR-29b, miR-15a, and a novel candidate whose target relationships persisted independently of the overall developmental trajectory. Notably, miR-29 showed particularly strong induction during maturation, consistent with its conserved roles in extracellular-matrix remodeling, differentiation, and vertebrate ageing. Together, these findings reveal coordinated remodeling of conserved miRNA-mRNA regulatory programs across the rapid maturation of N. furzeri and establish a framework for investigating post-transcriptional regulation during vertebrate development in a species with limited miRNA annotation.